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DEUTERIUM - HELIUM-3 FUSION POWER 
BALANCED CALCULATIONS 


by John J. Reinmann and Warren D. Rayle 
Lewis Research Center 

SUMMARY 

Steady-state operating parameters were calculated for the D-He cycle. Energy 
equations were written for the electrons and the two ion species in the plasma. The 
effects of varying amounts of reflection and reabsorption of cyclotron radiation were 
studied by the inclusion of a variable cyclotron radiation parameter. Other parameters 
were varied to assess the effects of (1) unequal confinement times of plasma electrons 
and ions, (2) plasma ion heating by fusion products at an enhanced rate, and (3) injection 
of fuel ions with varying amounts of initial energy. Results are presented in the form 
of equilibrium values of a containment parameter, nr, the electron and ion tempera- 
tures, and a breakdoAvn of the output power into bremsstrahlung, cyclotron radiation, 
and power carried out by escaping ions and electrons. Results are plotted as a function 
of plasma ion temperature. 


INTRODUCTION 

Power from controlled fusion reactors may be important in the future not only for 
ground but also for space applications (ref. 1). Present concepts of fusion propulsion 
systems, based on a steady-state D-He^ cycle, are presented in references 1 to 3. 

This cycle permits most of the fusion energy, carried by the charged reaction products 
to be retained within the magnetically confined plasma and also reduces the flux of ener- 
getic neutrons which may impose intolerable heat loads on superconductive magnet ele- 
ments. The energy carried by escaping charged particles, primarily the heated reac- 
tant ions, can be used directly to produce thrust (ref. 4). 

A better understanding of the probable range of operating conditions for such a 

O 

space -application reactor can be obtained by examining the energy balance of a D-He 
plasma. Energy balance studies of the DT cycle have been reported (see, e.g. , ref. 5) 



These equate the energy produced by fusion within the plasma to the net energy loss from 
the plasma and determine the conditions of density, confinement time, and ion and elec- 
tron temperatures which correspond to such an equilibrium. Studies envisioning ground 

applications emphasize the DT reaction, which has the lowest ignition temperature. The 
3 

D-He' reaction needed for space applications has not received equivalent attention. 

The purpose of the study reported here is to determine the required steady-state 

3 

operating parameters for the D-He reactor. The approach parallels that of refer- 
ence 5. Energy equations were written for the electrons and the two ion species in the 
plasma. Energy exchange rates were based on assumed Maxwellian distributions for 
these three species, as were the rates at which energy was added to each from the fusion 
product ions. Bremsstrahlung and cyclotron radiation losses from the electrons were 
also based on an assumed Maxwellian distribution, even through as Rose has pointed out 
there is reason to expect departures from such a distribution both at the very low and 
very high velocities. The cyclotron radiation parameter , C 2 /D, of reference 5, was 
varied through a range of values to represent the effects of varying amounts of reflection 
and reabsorption of the radiation within the plasma. Other parameters were varied to 
permit an assessment of the effects to be expected when (1) confinement times of elec- 
trons and ions are unequal, (2) ion heathy by the fusion products is at an enhanced rate, 

and (3) the reactant ions are injected with varying amounts of initial energy. 

4 

The effect of He and proton ash buildup was not considered in this initial study . 

Such ashes may be considered as having either of two effects: First, the positively 
charged He*^ and proton ashes will be accompanied by an equal number of electron 
charges. Hence, for a given fuel -ion density and temperature there will be proportion- 
ately larger bremsstrahlung and cyclotron radiation losses. Second, if the reactor den- 
sity is limited by the plasma pressure, the attainable density of fuel ions must be de- 
creased and the fusion power reduced. If the ash pressure is not too large , the results 

of this study should still be substantiaUy correct. 

3 

For this initial study, only the D-He reaction was included in the calculations. For 

ion temperatures greater than 30 keV, which applies to all cases in this report with zero 

3 

ion injection energy, the D-D reaction energy was less than 5 percent of the D-He reac- 
tion energy . Below 20 keV, the D-D reaction might become significant. For the case of 
energetic ion injection, steady -state operation at temperatures below 20 keV is possible . 
However , in most of these low temperature regions fusion power is insignificant com- 
pared to injection power . Therefore the trends predicted for energetic injection will 
still be valid. 

The results of the analyses are presented in the form of values of the containment 
parameter nr and electron and ion temperatures for which the energy balance is satis- 
fied over a range of values of the parameters. 
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PLASMA POWER BALANCE 


The work in this section is essentially the same as the power balance calciilations 
given in reference 5. The equations for fusion reaction rates were taken from refer- 
ence 6. The fast-ion energy-loss equations, and the plasma-component energy- 
exchange equations are in the forms given in the review article of reference 7 . The 
bremsstrahlimg and cyclotron radiation loss equations are from reference 5. Units are 
SI throughout, except as noted. 

The block diagram for the plasma power balance is given in figure 1 . The plasma 
components are electrons, deuterium ions and helium-3 ions, which are injected with 
energy Ej^, Ejj^, and Ejjj respectively. In all of the calculations of this report, Ej^ 
is set equal to zero. The energy distribution of all plasma components as assumed to 
be Maxwellian. The three plasma components exchange energy with each other, and 
they are heated by the fusion protons and « -particles - all by means of Coulomb colli- 
sions. The electrons lose energy by bremsstrahlimg Eg and cyclotron radiation E^. 
As plasma particles escape from the reaction chamber they carry their kinetic energy 
with them. 

O 

We assume a 50 percent D - 50 percent He plasma in the reaction chamber . The 
number density of deuterium ions, helium-3 ions and electrons is, respectively. 


nn = 


n 


( 2 ) 



(3) 


Since helium-3 has a double positive charge, there must be more electrons than 
ions in the plasma to maintain charge neutrality. The number density of protons and 
cv-particles is assumed to be negligible. 

The mean confinement times of the deuterium ions, helium-3 ions, and electrons 
are, respectively, Tg, t g, and t^. The proton and a -particle confinement times are 
assumed to be equal to their slowing down times as will be explained later . The escape 
rate of the s-component is WgAg- The mean kinetic energy carried by each particle is 

assumed to be (3/2) kT . 

s 
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Thermonuclear Reaction Rates 


Reaction rates (crv) for the deuterium - helium-3 reaction were taken from refer- 
ence 6. The empirical formula for (ov) , which Kozlov estimates to be accurate to with- 
in 15 percent, is 


(ov) = 0.90x10"^^ 3 

/l + 4 

where 

93.82 

and where T is in units of keV. 

In our calculations T was taken as (Tp + Tjj)/2. Since the deuterium and helium-3 
ions are strongly coupled, their temperatures are approximately equal for all the results 
in this report. 



Energy Loss from Fusion Protons and a-Particles 


The deuterium - helium-3 reaction is 

D + He^ — Proton (14.7 MeV) + He^ (3.6 MeV) 


The proton and cy-particle (He^ are born with energy far in excess of the plasma thermal 
energy, and for this reason they are referred to as fast ions. Before leaving the reaction 
chamber, the fusion particles will give up most of their energy to the three plasma com- 
ponents. The mean rate of change of kinetic energy of a fast ion k moving throv^h a 
plasma with N-components is (see, e.g. , ref. 7) 



s s s 




Vn ^ ' 



(5) 
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where s refers to the plasma components. 

All quantities are defined in the symbol list. The quantity h which does not appear 
in the usual slowing-down equation, is the enhancement coefficient. It is included to take 
into account the possibility that energy transfer from the fast ions to the plasma com- 
ponents may occur more rapidly than predicted by binary Coulomb collision theory . In 
all calculations, the Coulomb logarithm was set equal to 20. 

The energy transfer rate from the fast ion k to the plasma component s is defined 
as 


dE 


ks 


e^zj hm Z^L„ 
k s s s s 


dt 


^ks 


47re2vk 


m. 


$(XJ -_2_(l+rg)Xgexp(-x2 




( 6 ) 


Then let Rj^ be defined as 

N 

s=l 


{7) 


The ratio is the fraction of energy loss from the k ion that is absorbed by 

the s plasma component as the k ion energy is changed from to Ej^. The energy 

transfer from the fast ion k to the plasma component s is 



( 8 ) 


To determine the total energy transferred from the fast ion k to the plasma com- 
ponent s, we must select the upper limit of integration Ej^. As the energy of the fast ion 
decreases, it will reach a critical energy where there is zero mean energy ex- 

change between the k ion and the s plasma component, that is. 


dE 


ks 


dt 


= ■'ks = ° 


The critical energy 

s 

of the s plasma component 


depends on the mass ratio m^/mj^ and on the temperature 
Tg, Since there are three plasma components, there are 
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three 


values of . The largest of the three values of (Ej^^.^.) is used as the upper 


Kr) 


limit of integration for equation (8) . When the fast ion energy reaches 

/ \ / ^ “-'largest 

it is assumed to disappear from the plasma. Values of IT^ are taken from ref- 

erence 7 and appear in table I, Equation (8) is integrated numerically. 


TABLE I. - VALUES OF (%cr^s/'^s 


Plasma component 

Proton 

O' -particle 

Electron 

1.5 

1.5 

Deuteron 

.757 

1.17 

Helium-3 

.626 

1.065 


The fast ion slowing down time is defined as the time it takes to reach the energy 


E 


kcr 


. To obtain t. , we begin by combining equations (5) and (7) to give 


largest 


dt = 


dEj 


(9) 


After multiplying equation (9) by n, it can be integrated for nTj^ to give 



( 10 ) 


After the power balance equations are solved for the equilibrium values of Tj^, 
and Tg, the quantity nTj^ is calculated from equation (10) . The power balance equa- 
tions embody the assumption that the fast ion is contained within the plasma for a time 
eqiual to or greater than this slowing -down time Tj^. K the fast-ion confinement time is 
about the same as that of the plasma ions t the assumption can be tested by examining 

o 
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the ratio iiT„/nT, . The larger this ratio, the greater is the probability that the assump- 
tion is justified. The actual confinement time for the fast particles will depend on factors 
not included in the analysis. In some cases the confinement time may increase for fast 
ions because of their slower diffusion; in other cases the larger gyroradio may lead to 
more rapid loss rates. The ratio was calculated for the cases reported; when it 

drops below unity that fact is noted. 


Energy Exchange Between Plasma Components 


For the plasma components characterized by Maxwellian distributions, the energy 
transfer rate per unit volume from component a to component b is given as (see, e.g. , 
ref. 5) 



^('"a - ^b) 
’’e(a-b) 


keV 

(sec)(m^) 


( 11 ); 


where 


^e(a^b) - 


3(1.6021)^^^x10~^Q 

c^ 


^a”^b 

n^Lzhy 

b a b 



(12) 


where T^^ and T^ are in keV. The Coulomb logarithm L was set equal to 20 in all 
calculations. 

As pointed out in reference 5, the assumption of Maxwellian distribution for the 
electron population is not strictly sound. In particular, the energy transfer from the 
fast ions is primarily to those electrons with velocities less than the ion velocity. Con- 
sequently, the slow electrons are heated, and the distribution at low velocities depleted. 
This results in fewer electrons capable of receiving the energy from the fast ions, and a 
reduced energy transfer rate. Although this effect should be taken into account to obtain 
the best approximation to the energy balance, it is quite difficult to do so in a parametric 
analysis such as the present one. Moreover, the effect appears to be small at conditions 
of interest. The low velocity electrons diffuse rapidly in velocity space and tend to fill 
in any deficiency from the Maxwellian distribution. The assumption of Maxwellian dis- 
tribution for the electrons is also risky at the other extreme - the high velocity electron 
population will be depleted both by bremsstrahlung and cyclotron radiation losses as 
will be discussed later. 
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Bremsstrahlung Loss 


The bremsstrahlvu^ loss term, for a Maxwellian electron distribution, is 


E 


B 


3X10 n^ + njjZjj 




keV 

(sec)(m^) 


( 13 ) 


where isinkeV. 

In this initial study, normal bremsstrahlung was assumed, which corresponds to 
setting = 1 in reference 5. 

Cyclotron Radiation Loss 

The cyclotron radiative loss is the same as Rose's. It is 


Ec = C2X0.387 + 


keV 
(sec)(m’^)l 


where is the plasma transparency coefficient. 




2.1X10 T 


-3 rp 7/4 


if 


1/2 


The quantity if introduced in reference 5 is 


if, 

2 2 .3kT^ 


and for the D-He^ 


reactor case if 


becomes 


= lx 10^ 
2 


ZB/3 



(14) 


(15) 


( 16 ) 
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Substituting equations (15) and (16) into equation (14), we get the cyclotron radiative loss 
term 


E 


c ■ 








keV 

(sec)(m^)_ 


(17) 


where T is in keV. 

In the power balance equations that follow, a coefficient is used which is de- 
fined as 



(18) 


where I is the plasma size for an equivalent slab, ’’perhaps equal to the radius for 
cylindrical geometry” (ref. 5), B is the local magnetic field strength in the plasma, and 
j3 is the ratio of plasma pressure to the local magnetic field pressure. 

If the external magnetic field pressure balances the sum of the internal magnetic 
field pressure and plasma pressure, then can be written as 



where Bq is the external magnetic field strength, and /3 q is the ratio of plasma pres- 
sure to external magnetic field pressure. 

In reference 5, the quantity C 2 accounts for the fact that the cyclotron radiation is 

partly reflected at the electrically conducting vacuum wall. In reference 10, it was 

pointed out that the effect of reflectors is to increase the photon path length by a factor 

(1 - r)“^ where F is the reflectivity. Since photon path length is proportional to the 

system dimension I it is reasonable to assume (as is done in ref. 10) that the radiation 

loss is reduced by about (1 - F)^^^. Thus the quantity of reference 5 may be approx- 

1/2 ^ 

imated by the quantity (1 - F) ' . The preceding equation then becomes 
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(1 - 



Solving this last equation for F, gives 


r = 1 



(19) 


Figure 2 , which is a plot of the cyclotron reflectivity F against |3 q with as a 
parameter, is obtained from equation (19). The acceptable range of values of are 
obtained from energy balance considerations and from space propulsion requirements. 
Figure 2 will be referred to again in the RESULTS section. 

We have used approximation for as given in reference 5 without checking its 
validity in the range of electron temperatures and if values used in this study. Hope- 
fully, any errors should be adequately covered by the variable coefficient C 2 , which 
combines the effects of partial reflection of the cyclotron radiation at the vacuum cham- 
ber wall and the effects of K^. Departures from the assumed Maxwellian distribution 
may have particularly important effects on cyclotron radiation. In the high-velocity tail 
of the distribution the radiative losses are the greatest while the thermalization proceeds 
most slowly. The exact equilibrium will be determined by the balancing of the depopu- 
lating effects of cyclotron radiation and bremsstrahlung against the repopulating effects of 
diffusion in velocity space. (And even that balance should be modified to account for the 
probable velocity -dependence of the particle loss mechanisms.) The problem is made 
even more complex by the self -absorption of the cyclotron radiation. Energy is trans- 
ferred within the electron population as the electrons emit and absorb the (doppler- 
broadened) cyclotron fundamental and its many harmonics. Exact solutions for such a 
system seem to be beyond our present capabilities. One technique that may be used is to 
assume a net energy loss, for a given electron, equal to a fixed fraction of its single 
particle radiation rate. 
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Thermal Power Balance Equations 


With the equations presented in the previous sections it is a simple matter to write 
power balance equations for the three plasma components D, He , and electrons. First 
we write the particle balance equations for the deuterons and helium-3 ions, and obtain 
the expression for fuel fractional burnup. A word-form equation for the steady-state 
particle balance is 


Rate of change 


Injection 


Nuclear 


Geometrical 

of particles in 
a given volume 

— 

rate 


reaction 
loss rate 


escape rate 


The deuteron particle balance equation is 


dn, 


D 

dt 


Sjp - njjnjj(ov) 


n 


D 


'D 


= 0 


The helium-3 particle balance equation is 


( 20 ) 


( 21 ) 


dnrr nri 

— = Sjj -nj3njj<ov> 0 

at Tjj 


The fuel fractional burnup fg is defined as 


( 22 ) 


. _ Number of fuel particles reacting per unit time 



Number of fuel particles injected per unit time 



2 ngng<av) 
Sp, + Sr 


(23) 


Using equations (21) and (22) for and Sjj, respectively (noting that Tg = was 
assumed in all of the calculations) results in the following expression for fgi 



nTg(crv) 


(24) 
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Deuteron energy equation . - The steady-state deuteron energy equation is repre- 
sented in word-form as 


Net D 


D injec- 


a -fusion 


Proton 


Power 


Power 


Power 

heat- 


tion 


power 


fusion 


lost 


from 


from 

ing 

rate 


power 


toD 


power 
to D 

■ 

by escap- 
ing D 

- 

D to 
elec- 


D to 
He^ 

L J 










trons 




(25) 

Noting that the energy content of deuterons per unit volume is (3/2)npTj^ the deuteron 
energy equation becomes 

|(f "D^d) ^ Kt> * V - ^ J '^D 


2 '^e(D-e) 2 '^e(D-H) 


keV 

(sec)(m^). 


(26) 


where T and E are in keV. 

Helium-3 energy equation . - The steady-state energy equation for helium-3 is the 
same as equation (25) with the subscripts D and H interchanged. Therefore the 
helium~3 energy equation is 


d 3 


dt\2 


-hTh 




n 


H 


'H, 


‘H 3 


Eih + nj^njj(crv> (E^jj + E^jj) - ^ i Tjj 


Tr 2 


2 '^e(N^e) ^ '^e(N-D) 


keV 

(sec)(m^)_ 


(27) 


where T and E are in keV. 

Electron energy eqxiation . - In all calculations the electrons were assumed to be in- 
jected with negligible energy (Ej^ = 0). In word-form the energy balance equation for the 
plasma electrons is 
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Net 


a -fusion 


Proton 


Power 


Power 


Power 


Brems- 


Cyclo- 

elec- 


power to 


fusion 


lost by 


from 


from 


strahlung 


tron 

tron 


electrons 


power 


escap- 


elec- 


elec- 


power 


power 

heat^ 




to elec- 


ing 


trons 


trons 


loss 


loss 

ing 




trons 


elec- 


to D 


to He^ 





rate 






trons 










The electron energy equation is 


dt 



= nj^njj(av> 


(E + 
' ae 


Epe) - 


3 npnjj(crv) 


+ 3 -j. _ 3 ^e^'^e ~ 

^e/ ^ ^ ^ ^e(e-D) 


3 ^e^'^e ' 


- Eb - Ec = 0 


keV 

(sec)(m^)_ 


where T and E are in keV. 


(29) 


Solution of Three Simultaneous Energy Equations 

We make use of equations (1) , (2) , and (3) to get a common number density n in 
the three energy equations. Equations (26), (27), and (29) can be written in the following 
functional forms; 


— =«Td. t,. Ejo) 


nr- 


D 


(26a) 


nr 


/’’h’ '^e’ 


(27a) 


H 


nr^ 


h(Ti,, Tg, T^, cy 


(29a) 


There are also the assumed relations 
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iitd = nTjj 


(30) 


and 


nxg = Jnxj^ (31) 

Equations (26), (27), (29), (30), and (31) are five simultaneous equations in terms of the 
10 variables nr^, nTjj, nr^, ’^D’ "^e’ ^ID’ ^IH’ ^2’ and J. Therefore, the quan- 

tities Ejjj, Ejjj, C^, J, and are chosen as input parameters, and the variables to 
be determined are nvj^, nxjj, nr^, T^, and Tjj. These five simvdtaneous equations 
can be reduced to two simultaneous equations with T^ and Tjj as unknown variables. 

Since, from equation (30), = Tjj, we can eliminate nvj^ from equations (26) 

and (28). After some rearrangement the following equation results: 


( OV*) (^jjj 


+ E TT 

q;H 


Epn) 


3(Th - T,) 
”'^e(H-*e) 


3(Th - Tp) 

’^'^e(H-D) 


3 Tj) - 2 Ejjj 


X 


<crv>(Ejj5 + E^j3+ Epj^) - 


3(Tp - T^) 
'''■e(D-e) 


- T^) 
^'^e(D-H)_ 


(32) 


Use equation (31) to replace nr^ by JnVj^ in the electron energy equation (29). 
Then eliminate nvj^ between equations (26) and the electron energy equation (29) to ob- 
tain the following result. 




9(T^-Tp) 9(T^-Tg) 4 Eb 4Eg 9 

’^'^e(e-D) ^■^e(e-H) n^ n^ ^ 


X 


(^>(E^p + Epj5 


Eid) 


3(Tp - T^) 
^^e(D-.e) 


3(Tp - Th) 
^^e(D-H) 


= 0 


(33) 


Equations (32) and (33) are two simultaneous equations with Tjj and Tjj as un- 
knowns, and Ejj^, Ejjj, C^, J, and T^ are inputs. The solutions, Tj^ and Tjj, are 
substituted into equation (26) to obtain the corresponding value of nxjj. Similarly, the 
values of fast-ion slowing down times nvj^, fractional fuel burnup f^^, and the various 
forms of output power are all calculated after Tj^ and Tjj are obtained. When Ej^ 
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and Ejjj are zero, there is only one solution for Tj^ and Tjj. When Ejp and Ejjj 
are not zero there are sometimes two solutions. Solutions are obtained on a high-speed 
digital computer. 


RESULTS 

Zero Ion Injection Energy 

The cases with zero injection energy have the results shown on figures 3 to 5. Fig- 
ures 3(a), 4(a), and 5(a) are for the normal or most probable case, with equal confine- 
ment times for electrons and ions (t^ = Tj^) and with heating of the ions and electrons by 
the fusion product ions occurring at the basic rates calculated from binary collisions 
(he = hj^ = hjj = 1) . Figures 3(b) , 4(b) , and 5(b) explore the effect of assumed enhanced 
fast-ion heating of the deuterons and helium-3 ions. The deuterons and helium-3 ions 
are assumed to gain energy from the fast ions at 1000 times the normal rate (h^ = hjj = 
1000, hg = 1). Finally, figures 3(c), 4(c), and 5(c) return to normal ion heating but as- 
sume the electron confinement time to be 10 times that of the ions. The opposite as- 
sumption, that the confinement time of the electrons is one-tenth that of the ions, gives 
the results shown on figures 3(d), 4(d), and 5(d). Throughout f^ures 3 to 5, the results 
are plotted as a function of deuteron temperature. The cyclotron radiation coeffi- 
cient C 2 varies as a parameter from zero to 2. 5. 

Two points to keep in mind throughout this discussion are that (1) the confinement 
times of the deuterons and helium-3 ions were assumed to be equal (i.e. , = Tjj) and 

(2) because of the strong coupling between the deuterons and helium-3 ions, their tem- 
peratxires were always approximately the same to within 10 percent or less. 


Confinement Parameter, nxp 

Figure 3 shows the confinement parameter nxj-j as a function of deuterium ion tem- 
perature. This relation alone does not determine the operating conditions for any par- 
ticular reactor . The reactor will have a confinement geometry which will impose some 
relation between n, Tj^, and T. If this relation can be expressed as nr^j = f(T), to take 
a very simple example, the actual operating points would be represented by the intersec- 
tions of such a curve with the energy balance nr curve. In some respects, such a sys- 
tem can be likened to the case of the electronic vacuum tube, where the tube characteris- 
tics must be combined with the load line representing the external circuitry. 
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There are several interesting pieces of information on figure 3(a). First, the min- 

20 

imum value of is 7x10 seconds per cubic meter, for zero cyclotron radiation 

loss. This is about a factor of five above the minimum value for a closed DT system. 
Second, the minimum ion temperature for steady-state operation is about 30 keV. Near 
this point, the particle loss rate approaches zero. For lower ion temperatures, the 
bremsstrahlung loss alone exceeds the fusion power released, precluding steady-state 
operation (without energetic ion injection). Third, the effect of increased cyclotron ra- 
diation loss is to increase the minimum nr^ and to impose an upper limit on ion tem- 
perature. This upper limit represents the case wherein the total energy produced by 
fusion is lost by the bremsstrahlung plus cyclotron radiation, leaving none to escape with 
particles. At higher temperatures, the cyclotron radiation increases more rapidly than 
the fusion power so that the radiative losses exceed the fusion power. Hence steady- 
state operation at higher temperatures is not possible. 

All four cases shown in figure 3 exhibit the same trends as the cyclotron radiation 
parameter, C^, is varied. Figure 3(b) shows the effect of enhanced heatii^ of the plas- 
ma ions by the fusion products. Since less energy goes directly to the electrons, the 
electron temperature should be suppressed and the radiative losses reduced particularly 
at the higher ion temperatures. This effect is indeed noticed. The upper limit on ion 
temperature is increased over the normal case, pronouncedly so at the higher values of 
C^- The minimum values of nr are lower than the normal case: only slighter lower 
for C 2 = 0 but lower by a factor of 10 for = 2. 5. Any enhanced ion heating effect 
would be very welcome. 

Although it seems likely that in most reactor schemes the confinement times of the 
electrons and ions will be about equal, occasionally the idea of unequal throughputs of 
electrons and ions is considered. Figure 3(c) shows results obtained when = 10. 

The values of nxjj are slightly less than the normal case, because the energy lost with 
the escaping electrons is less, which must be compensated by increased ion loss rates. 
However, the changes are small. Even if one could devise a means for preferentially 
confining the electrons, the incentive seems inadequate. 

Figure 3(d) represents the opposite case = It seems more plausible, 

from an engineering viewpoint, to consider increasing the flow of cold electrons into the 
plasma by a factor of 10. The first thing we note is that nxjQ for the ions is about a 
factor of five greater than for our normal case, while the allowable range of ion temper- 
ature is about the same . The effect may be viewed as an increase in required ion- 
confinement time; hence, a reduction in energy transported out with the ions which com- 
pensates for the increased energy lost with the increased electron flow. 
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Electron Temperature 


The values of electron temperature calculated for the four cases are shown on fig- 
ure 4. The normal case (fig. 4(a)) shows that, for zero cyclotron radiation, the electron 
temperature slightly exceeds that of the ions up to about 175 keV. At still higher tem- 
peratures, the combination of increased bremsstrahlung losses and proportionately less 
fast-ion heating of the electrons causes the electron temperature to fall below that of the 
ions. Increased cyclotron radiation losses act to cool the electrons, especially at the 
higher temperatures. 

In figure 4(b) , the effect of enhanced ion heating is to transfer proportionately more 
fusion energy directly to the plasma ions so that the electron temperature is seen to fall 
below that of the ions. Most of the electron heating now comes from the plasma ions, 
rather than from the fast fusion products. This was confirmed by calculations setting 
hg = 0, that is, no electron heating directly from the fast ions. The resulting electron 
temperatures, for the = 0 case, were approximately the same as in figure 4(b). 
Energy exchange with the plasma ions is sufficient to keep the electron temperature from 
dropping below about one-half of the plasma ion temperature. 

Figure 4(c) demonstrates the effect of increased confinement of the electrons. For ; 
zero cyclotron radiation, the electron temperature rises sharply above the ion tempera- 
ture. This permits the increased loss by bremsstrahlung and the increased energy lost 
with each escaping electron to partly compensate for the reduced number of electrons 
escaping. Because of the low values of nvj^, for the case of Cg = 0, the ratio of the 
plasma ion confinement time to the fast-proton confinement time is less than unity for 
Tj^ above 62 keV. However, the ratio is always greater than one -third, so that the elec- 
trons at least are confined longer than the proton slowing down time. The results for 
Cg 0 are much closer to the corresponding normal cases. The cyclotron radiation 
losses are quite sensitive to electron temperature, so that a small increase in electron 
temperature serves to restore the balance at the higher C^- 

Results for the case of rapid electron throughflow are plotted on figure 4(d) . With 
Tg/Tp = 0.1, the electron temperature is strongly suppressed. The large flux of elec- 
trons serves to carry off most of the energy the electrons receive, even at a lower value 
of the mean energy per electron. 


Distribution of Output Power 

In the introduction of this report it was noted that for a fusion propulsion system, 
the majority of fusion energy must appear in the escaping plasma particles. When these 
particles are emitted preferentially in one direction, they produce thrust. This thrust 
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may be augmented by using the escaping particles to ionize and heat additional propel- 
lant, followed by an expansion in a magnetic nozzle. A theoretical analysis of such a 
process (ref. 4) revealed that the plasma ions of our energy rar^e could transfer energy 
effectively to the propellant, but that electrons of equal energy were relatively ineffec- 
tive. We are interested, then, in how the total energy output from a fusion reactor is 
distributed among the various forms, that is, in bremsstrahlung and cyclotron radiation, 
in escaping electrons and in escaping ions. These quantities have been calculated and 
are presented in figure 5. The plots in figure 5 are the fraction of fusion energy re- 
leased to the plasma that appear in various forms of output energy. 

The results for the normal case, with zero cyclotron radiation, are presented in 
figure 5(a-l) . Although the ideal case of zero cyclotron radiation is unattainable, it es- 
tablishes an upper limit on the fraction of the output in the form of charged particle kine- 
tic energy. At best about 84 percent of the output will be in the form of charged parti- 
cles. The electron and ion temperatures for this case are approximately equal, so that 
the plasma components carry out energy in proportion to their number densities. The 
electrons carry about 1. 5 times the energy of the ions. Also, as Td is reduced to 
about 30 keV, all of the energy appears as bremsstrahlung. At the higher temperatures, 
the bremsstrahlung fraction is nearly constant. The bremsstrahlung increase with 
Tg^^^ is nearly matched by a corresponding slow increase in the fusion reaction rate 
(ctv) . The effect of cyclotron radiative loss on the distribution of output power is shown 
in figures 5(a-2) to 5(a-4) for between 0. 1 and 2.5. The fraction of power output in 
charged particle kinetic energy varies with ion temperature, with the maximum fraction 
shifting to lower ion temperatures as is increased. The maximum assumed cyclo- 
tron radiation loss (C^ = 2.5) would result in very poor operation, for only about 11 per- 
cent of the energy output is in charged particles. Worse, only about 4. 4 percent is in 
the ions. 

We conclude two things from the results thus far. First, large values of cyclotron 

radiation imply higher confinement times; second, for the normal case of equal electron 

and ion confinement times the electrons carry away a disproportionate share of the fusion 

energy. Perhaps, as is common in open-ended systems, some of this electron kinetic 

energy can be transferred to the ions by means of an ambipolar potential drop at the 

exit of the reactor. Third, if at least 40 percent of the output energy is to be in the form 

3 

of charged particles, must be less than 1.0, and the D-He reactor must operate 
somewhere between 40 to 70 keV. 

In our attempts to increase the fraction of energy carried out by the plasma ions, 
we considered two extreme variations from the normal conditions. First, it may be pos- 
sible to get enhanced heating of the plasma ions by the fast ions. Figure 5(b) shows the 
result of assuming 1000 times the normal energy transfer rates between these species. 
(Results are substantially the same for any enhancement factor of 1000 or greater.) 
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Figure 5(b-l) is again our unachievable ideal (C^ = 0). The electron temperatures are 

lower, thus the bremsstrahlung fraction is now only 11 percent at the best condition. 

For the same reason, the energy carried by the electrons has dropped. Now the ions 

carry about 1.25 times as much as the electrons. The effect of increasing Cg is much 

less pronounced than in the normal case because of the lower electron temperatures. 

These results are very attractive, even for the worst case of Ck = 2. 5, the charged 
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particles carry about 47 percent of the energy. The minimum uTj^ is about 9x10 
seconds per cubic meter. At present, we know of no physical process which would pro- 
vide enhancement of this magnitude. 

Another scheme for getting a higher fraction of the fusion output into the escaping 
ions is by reducing the electron loss-rate relative to the ion loss rate. This scheme 
also is studied in figure 5(c) for the case of = 10. Again, for the ideal = 0 

(fig. 5(c-l)) we see that indeed the distribution of energies is good. Although the frac- 
tion of bremsstrahlung has increased to about 22 percent minimum, the plasma ions es- 
cape with about 1.4 times the energy of the escaping electrons. Unfortunately, increas- 
ing Tg also causes an increased electron temperature. Hence, the radiation losses are 
large, causing the fraction of power output in escaping particles to be low at high 
values. If values of on the order of 0. 1 are achievable, then there would be some 
advantage to this scheme, namely, more energy output in the ions. 

For completeness, we include the results of the case where the electron throughput 
is 10 times the ion throughput in figure 5(d). One might hope for some encouraging re- 
sults because the electron temperatures are suppressed. However, looking at the 
C 2 = 0 case, figure 5(d-l) we see that althov^h bremsstrahlung is slightly lower than the 
normal case, the fraction of energy carried out by the escaping electrons is about 
10 times that carried by the ions. This trend persists as Cg is increased. Since cy- 
cloton radiation losses are still very high as is increased, there seems to be nothing 
to recommend increasing the ratio of ion confinement time to electron confinement time. 


Reactor Beta Requirements 

It was pointed out in a previous section that must be less than 1. 0, for the nor- 
mal case, while for enhanced fast-ion heating, values of up to 2. 5 might be accep- 
table. Some insight into the consequences of these requirements on can be obtained 
from figure 2. Figure 2 is a plot of the reflectivity for cyclotron radiation F against 
|3 q with C 2 as a parameter. The product ZBq is taken to be 10.0. For the normal 
case where must be less than 1.0, values of F must lie above the curve for 
C 2 = 1.0. Hence, at very low values of |3 q near-perfect reflectivity is required. For 
/3 q greater than 0. 1, less stringent values of F are required. Values of ZBq greater 
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than 10.0 wotdd also help. Actual values of T for cyclotron radiation are not presently 
available, but a possible method for determining T has been presented in reference 8. 
Clearly low operation demands the difficult requirement of obtaining a surface, op- 
erated at extremely high material temperatures, with a high reflectivity to cyclotron 
radiation. On the other hand, if r is low, there is the equally difficult requirement of 
operating a closed reactor at high /3 q (ref. 9). 


Energetic Ion Injection 

Next, we explore the effects of energetic ion injection on steady-state operation of a 

3 

D-He reactor. For space propulsion application, it seems doubtful that energetic ion 
injection would be acceptable. The energy conversion system and other appurtenances 
required would so increase the system mass as to make it noncompetitive with other ad- 
vanced propulsion systems. However, the results with injection may be of interest to 
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those considerii^ the D-He cycle in conjunction with direct energy conversion concepts. 

The results, with various assumed ion injection energies, are displayed on figure 6, 
where uTj^ is plotted against deuterium ion temperature . The electron temperature 
contours are also shown. These results were obtained assuming zero cyclotron radia- 
tion, normal bremsstrahlui^, equal electron and ion confinement times and cold electron 
injection. One of the interesting effects of energetic ion injection is that it significantly 
expands the range of possible steady-state operation. In fact, for some combinations of 
electron temperature and ion injection energy there are two different operating points 
possible: One is at a high value of nTjj and low ion temperature; the other occurs at a 
lower nTj^ and a hi^er ion temperature. 

Some of the features of this map can be readily understood in physical terms. It is 
helpful to think of the density n as being held constant so that variations in nr^ imply 
variations in only. Figure 6 can be divided for discussion into three regions. First 
consider the upper-middle and upper -right region. At the higher operating temperatures, 
the curves for ion injection energies from 0 to 90 keV are essentially alike. Fusion re- 
mains the dominant energy source, but the addition of energy via injection permits oper- 
ation with lower confinement times. Near 30 keV the nr^ again increases without 
bound; here the injection power and the charged particle output power become negli- 
gible, and a balance is found between the bremsstrahlung and the fusion terms. At the 
higher injection energies, as Tj^ decreases the curves turn down instead of up and pro- 
vide operating conditions at low values of nTjj which will be discussed later . 

In this upper -right region the ratio of ion confinement time to fast-ion slowir^ down 
time is often less than one because of the lower values of nr^j. Hence the numerical 
values of nTjQ in these regions are less reliable, but still are useful for indicating 
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trends. If less fusion energy is absorbed by the plasma than the calculations indicate, 
the nxj^ curves would be squeezed together and pushed upward toward the zero injection 
case. In those cases where the higher injection energy curves turn downward, the fusibh 
power becomes negligible compared to the injection power so that the^ curves should not 
be significantly changed by the fact that the slowing down time of the reaction products is 
very much greater than the ion confinement time . ■ 

In the upper-left part of figure 6 is a set of curves heretofore absent, but of little 
practical significance. For very low ion temperatures and very long confinement times, 
the injection energy term is low. The fusion term is much lower. Consequently, a 
small injection energy is balanced by a small bremsstrahlung with electron and ion tem- 
peratures equalized at a low value. 

In this same region, near 30 keV ion temperature the confinement times are long 
and injection power is small compared with fusion power. Again, the fusion energy is 
balanced against bremsstrahlung losses. , 

The last region to consider is the lower portion of figure 6. Here the throughput of 
ions and electrons becomes so large that the injection energy and the particle loss terms 
dominate. It m^ht correspond to an experimental machine with very poor confinementi 
As the injected ions pass through the machine, they become thermalized (according to 
our model) and exit with a temperature approaching 2/3 of the injection energy. These 
possible operating conditions have little practical importance, because the ratios of re- 
circulating power to fusion power become very large. For example, with a 90-ke¥ in- 
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jection energy, this ratio has already reached 10 at an nr^ of about ixlO seconds per 

cubic meter. For a 180-keV injection energy, the ratio reached 10 at about nxjj = 

2.7x10^^. These are high values of nr^ to expect in an open-ended machine. 

Other assumptions in our model start to break down as nxp, becomes low. For 
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values of nr^ less than 10 or 10 some cases will not permit confinement times 
long enough to thermalize the plasma ions, or to slow down the fast ions and transfer 
their energy to the plasma. Although the model could perhaps be modified to account for 
the resulting non-Maxwell ian distributions and partial energy transfer, it seems hardly 
worth the effort. 


CONCLUDING REMARKS 

3 

Solutions to the energy balance equations for the D-He fusion plasma have been ob- 
tained. These indicate possible operating conditions required for such a fusion reactionl 
The validity of these solutions for an actual plasma depends on the validity of the many I 
assumptions embodied in the original equations. Among these assumptions are (1) that 
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the electrons, deuterons, and He ions within the plasma are characterized by a 
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Maxwellian distribution, (2) that the energy interchanges among plasma species as well 
as the fast -ion heating could be derived from the simplest, binary interaction (3) that 
the net cyclotron radiation loss may be represented by the approximate expression of 
reference 5, and (4) that the confinement time of the fast ions produced by fusion will 
about equal that of the plasma ions, or at least will exceed the time required for their 
energy to be lost to the plasma. 

The results show a minimum value of the containment parameter nr of about 
7 k 10 seconds per cubic meter for the normal case with neither injection energy nor 
cyclotron radiation loss. This is about five times greater than the value for the DT 
reaction, and occurs at a high value of about 100 keV ion temperature. When cyclotron 
radiation loss is included, the minimum value of nr is greater still and solutions exist 
only within a band of ion temperatures. The lower limit of this band, about 30 keV, is 
relatively insensitive to the magnitude of the cyclotron radiation loss. However, if the 
coefficient exceeds a value of about 2. 5 steady-state operation is no longer possible 
at any ion temperature. 

The minimum required nr may be reduced if it is possible to produce enhanced 
fast-ion heating. The reduction is particularly dramatic when the cyclotron radiation 
coefficient Cg is large. Lower values of nr may also be obtained by resorting to 
energetic ion injection. Such solutions are of limited interest -when nr is reduced by 
more than about a factor of 10 the injection power exceeds the fusion power by a like 
factor. 

For space propulsion or for direct conversion applications, the fraction of the total 
power carried by the escaping charged particles is important. When cyclotron radiation 
is a factor, this fraction exhibits a maximum at a particular value of ion temperature. 
The magnitude of this maximum varies from 0.1 to 0.7 for the normal case over the 
range of cyclotron losses considered. To maximize the particle energy losses, we must 
select an operating temperature just slightly below that for which nTj^ is minimum. 

The optimum operating temperatures range from 65 to 45 keV as Cr, varies from 0. 1 
to 2.5. 

The big question remains that of determining the appropriate value for the cyclotron 
radiation losses. At the higher operating temperatures the electrons radiate very 
strongly and over a large number of harmonics. At the same time, the plasma density 
and the reabsorption rate is reduced if we hold the plasma pressure constant. On the 
other hand, at these high temperatures the loss of energy by radiation will tend to distort 
the electron distribution and reduce the number of energetic electrons which radiate. If, 
for example, we require that 40 percent of the fusion power be carried by the escaping 
charged particles, the cyclotron radiation coefficient Cg has to be 1. 0 or less. If we 
could obtain enhanced ion heating by a factor of 1000, we would meet this goal with 
as great as 2.5. 
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Unless a suitable cyclotron-radiation reflector is developed (r >0.9), the D-He^ 
reactor will have to operate with /3q in excess of 0.1. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, January 22, 1971, 

129-02. 
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APPENDIX - SYMBOLS 

magnetic field 

cyclotron radiation loss coefficient 

C2/D 

velocity of light 

bremsstrahlung power loss per unit volume 

cyclotron radiation power loss per unit volume 

injection energy of s plasma component 

energy transferred from fast ion k to plasma component s 

initial energy of k fusion reaction product 

energy of k fusion reaction product at time t 

energy of fast ion k at which the mean energy exchange rate with the 
s plasma component is zero 

electronic charge 

fractional fuel burnup 

fast-ion heating enhancement coefficient 

ratio of electron to ion confinement times 

plasma transparency coefficient for cyclotron radiation 

Boltzmann's constant 

Coulomb logarithm 

Coulomb logarithm for s plasma component 
cyclotron radiation parameter (see eq. (15)) 
plasma size for an equivalent slab 
mass of s plasma component 
mass of k fusion reaction product 
number density of s plasma component 
total ion density 
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energy transfer rate per unit volume from plasma component a to plasma 
component b 

Rj^ mean energy transfer rate per unit volume from fast ion k to all plasma 

components 

^ks mean energy transfer rate per vmit volume from fast ion k to plasma 

component s 

Sg particle injection rate per unit volume of s plasma component 

Tg temperature of s plasma component 

t time 

Vj^ velocity of k fast ion at energy 

X dimensionless parameter 

Zg charge number for s plasma component 

Zj^ charge number lor k fast ion 

|3 ratio of plasma pressure to magnetic pressure 

r surface reflectivity for cyclotron radiation 

■)/g mass ratio mg/mj^ 

permittivity of free space niunerically equal to 8.854x10“^^ coul^/n-m^ 

9 normalized plasma ion temperature (see eq. (4)) 

T mean particle confinement time ; 

Subscripts: 

D deuterium ion 

e electron 

H helium -3 ion 

I injection 

k fusion reaction fast ion 

p proton 

s plasma components 

a helium -4 ion 
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Ratio of plasma pressure to external magnetic field pressure, pQ 


Figure 2. - Required reflectivity, r, for cyclotron radiation as function of 
external beta, pQ. 
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(a -4) Cyclotron radiation coefficient, 2.5. 
(a) Concluded. 

Figures. - Continued. 
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Power fraction 




(b) Confinement time ratio Xg/XQ, 1; enhanced ion heating (1000 times normal). 
Figure 5. - Continued. 
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Deuterium ion temperature, Tq, keV 
(b-4) Cyclotron radiation coefficient, 2.5. 


(b) Concluded. 
Figures. - Continued. 
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r fraction 




(c~2) Cyclotron radiation coefficient, 0. 1. 

(c) Confinement time ratio tg/iQ,. 10; norma! ion heating. 
Figure 5. - Continued. 
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(d) Confinement time ratio tg/iQ, 0. 1; norma! ion heating. 
Figure 5. - Continued. 
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Deuterium ion temperature, Tp, keV 
(d-3) Cyclotron radiation coefficient, 1.0. 

(d) Concluded. 

Figure 5. - Concluded. 



Containment parameter, ni^, sec/m 


Ion injection 
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